广 上 VAN 人 于 甘 日 于 || 
UnInanAIV 吕 上 只 有 下 
第 39 卷 第 8 其 计算 机 应 用 研究 TVol 39 No. 8 
录用 定稿 Application Research of Computers Accepted Paper 


一 种 基于 脑 电 微 分 蚁 的 警觉 度 估计 方法 
陈 ， 万 ， 华 艳 乎 ， 李 爱 华 ， 杨 梅 枝 ， 姜 ” 柯 
(火箭 军工 程 大 学 ,305 教研 室 , 西安 710025) 


摘 要 : 为 提高 人 员 警 觉 度 实时 估计 精度 ， 提 出 了 一 种 基于 微分 灶 (DE)、 改 进 滑动 平均 和 双向 二 维 主 成 分 分 析 (TD- 
2DPCA) 的 警觉 度 实时 估计 方法 。 首 先 将 总 频带 以 茶 一 带宽 分 解 为 多 个 子 频 带 ， 在 每 个 子 频 带 提取 DE; 然后 结合 警 
觉 度 的 时 间 动 态 特征 对 传统 滑动 平均 方法 进行 了 改进 ， 利 用 改进 滑动 平均 对 DE 进行 了 平滑 处 理 ; 之 后 利用 TD- 
2DPCA 对 DE 进行 降 维 ， 并 采用 最 小 二 乘 支持 向 量 机 (LS-SVM) 建 立 特征 抵 阵 与 警觉 度 之 间 的 回归 模型 ， 以 实现 警 
觉 度 的 实时 准确 估计 。 最 后 利用 SEED-VIG 数据 集 进 行 了 实验 验证 ,结果 表明 , 采用 改进 滑动 平均 和 TD-2DPCA 方 
法 对 数据 进行 处 理 可 以 提高 警觉 度 估计 的 精度 并 降低 估计 时 间 ; 总 频带 在 0-35HZ 以 内 ， 分 解 带宽 为 1HZ 或 2HZ 
时 ,提取 的 DE 进行 警觉 度 估计 能 够 达到 最 高 的 估计 精度 ，Pearson 相关 系数 约 为 0.91，RMSE 约 为 0.09， 优 于 现 有 
的 警觉 度 估计 方法 。 
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Vigilance estimation based on differential entropy of EEG 


Chen Wan, Cai Yanping', Li Aihua, Yang Meizhi, Jiang Ke 
(Room 305, Rocket Force University of Engineering, Xi*an 710025, China) 


Abstract: To improve the real-time estimation accuracy of vigilance, this paper developed a vigilance estimation method 
based on differential entropy (DE) , improved moving average and two-directional two-dimensional principal component 
analysis (TD-2DPCA) . First, it decomposed the total frequency band into multiple sub-bands with a certain bandwidth, and 
extracted DE from the sub-band. Then combined with the temporal dynamic characteristics of vigilance, it improved the 
traditional moving average and used improved moving average to smooth the DE. After that, it used TD-2DPCA to reduce the 
dimensionality of DE, and established a regression model between the feature matrix and vigilance with least squares-support 
vector machine (LS-SVM) to achieve accurate real-time estimation of vigilance. Finally, this paper used the SEED-VIG for 
experimental verification. The results show that the improved moving average and TD-2DPCA can improve the accuracy of 
vigilance estimation and reduce the estimation time. Using 1HZ or 2HZ bandwidth to decompose the total frequency band 
within 0-35HZ to extract DE and perform vigilance estimation, the result has the highest accuracy. The Pearson correlation 
coefficient is about 0.91, and the RMSE is about 0.09, which is better than the existing vigilance estimation methods. 
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中 采用 (9+0)1B 的 功率 谱 密 度 比 值 、 专 注 度 与 放松 度 的 比值 


0 引言 和 膀 眼 频率 作为 疲劳 检测 特征 值 ， 实 现 了 人 员 疲劳 的 实时 检 

警觉 度 是 指 人 集中 注意 力 执行 革 项 任务 时 所 表现 出 的 灵 “， 测 ， 然 而 该 方法 只 将 疲劳 状态 分 为 三 类 ， 无 法 实现 人 员 警觉 
敏 程度 ， 警 觉 度 越 低 ， 人 员 在 执行 任务 的 过 程 中 犯错 的 概率 度 的 准确 监测 。 张 朋 等 加 在 每 个 通道 信号 上 提取 了 37 个 特征 
就 会 越 大 中 。 生 活 中 存在 大 量 危 险 性 较 高 的 职业 ， 比 如 长 途 。 ”利用 随机 森林 方法 实现 了 人 员 的 疲劳 检测 ， 然 而 该 方法 没有 
客车 驾驶 员 、 化 工 品 运输 司机 、 大 型 设备 操作 员 、 飞 行 员 等 ， “对 提取 的 特征 进行 处 理 , 导致 特征 维度 较 高 , 识别 精度 有 限 。 
都 需要 人 员 长 时 间 保持 高 敬 觉 度 状态 ， 一 旦 人 员 出 现 警觉 度 。 ”上述 两 种 方法 提取 的 脑 疲劳 特征 都 是 来 自 EEG 信号 的 五 个 
降低 ， 可 能 导致 灾难 性 后 果 Dl。 因此 , 进行 人 员 警 觉 度 的 估计 “频段 信号 , 无 法 最 大 化 EBG 包含 的 信息 。 而 采用 经 验 模 态 分 
已 经 成 为 人 机 交互 和 主动 安全 领域 的 重要 研究 课题 。 解 、 小 波 变 换 等 分 解 方法 对 EEG 信号 进行 分 解 mg, 可 以 提 

常用 于 人 员 警 觉 度 佑 计 的 方法 有 自我 估计 、 眼 球 追踪 测 。 “取出 更 多 有 用 的 特征 ， 但 信号 分 量 含义 不 明显 ， 同 时 多 余 的 


量 、 脑 电 图 (Electroencephalogram, EEG) 测 量 , 其 中 基于 EEG 特征 也 可 能 对 原 有 特征 带 来 和 干扰， 影响 估 计 结 果 00。 同 时 ， 
的 人 员 警 惕 度 估计 方法 表现 出 更 高 的 准确 率 和 稳定 性 中 。 人 的 疲劳 状态 变化 是 一 个 动态 过 程 ， 不 能 简单 地 将 人 的 警觉 
EEG 信和 号 是 记录 头皮 电位 变化 的 信号 , 在 一 定 程度 上 反映 了 度 划 分 为 几 类 ， 因 此 实现 人 员 和 警觉 度 的 实时 评估 ， 更 加 贴 合 
人 的 认 知 活动 由。 通常 将 EEG 信号 划分 为 五 个 频段 ， 分 别 是 实际 应 用 [9 。 研 究 表明 03449, 采用 微分 粒 (differential entropy， 
0(0.5-4HZ)、0(4-8HZ)、2% (8-14HZ)、B(14-30HZ)、7 (30- DE) 作 为 警觉 度 估计 特征 效果 优 于 传统 的 脑 电 特征 ， 因 此 本 
50HZ)， 不 同 实验 中 划分 略 有 不 同 喇 。 研 究 表明 ， 昏 睡 状 态 下 文采 用 DE 作为 警觉 度 估计 特征 。 
0 和 wx 波 的 能 量 会 增加 ，4 和 7 波 的 能 量 会 降低 名 。 王 和 任 等 在 现 有 研究 基础 上 ， 本 文 进 一 步 对 人 员 警 觉 度 估计 方法 
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进行 了 研究 。 在 特征 提取 上 , 先 将 每 个 通道 EEG 信号 等 分 为 从 zavnn 开始 由 前 往 后 计算 平滑 结果 : 
多 个 子 频 带 , 然后 在 每 一 个 子 频带 上 提取 DE, 并 分 析 了 不 同 1 如 
等 分 方法 对 警觉 度 估计 的 影响 ， 以 期 找到 最 优 的 DE 提取 方 a 
法 ， 同 时 与 传统 特征 提取 方法 进行 了 对 比 。 在 特征 平滑 上 ， i 
结合 警觉 度 的 时 间 动 态 特性 ， 对 传统 的 滑动 平均 方法 进行 了 2 Dt pn 
改进 ， 提 出 了 变 窗 口 平 滑 与 短 时 预测 相 结 合 的 边缘 数据 平滑 第 三 步 : 由 于 区 间 [s] 和 区 -s+1 刀 对 应 的 窗口 太 小 ， 求 
方法 ,提高 了 边缘 数据 与 整体 数据 的 关联 性 。 在 数据 降 维 上 ， 平均 后 的 数据 仍 不 够 平滑 ， 因 此 采用 平滑 后 的 数据 进行 短 时 
对 比分 析 了 传统 主 成 分 分 析 (Principal component analysis ， 预测 计算 : 
PCA) 方 法 和 双向 二 维 主 成 分 分 析 方 法 在 脑 电 特征 降 维 中 的 = j=1,2,...,5 (5) 
性 能 。 最 终 ， 利 用 最 小 二 乘 支持 向 量 机 (Least squares-support zt(za -za)/n j=L-stl,...,L 
vector machine，LS-SVMD) 构 建 了 特征 向 量 与 警觉 度 之 间 的 匠 其 中 : :和 是 两 个 变量 ， 需 要 根据 实际 处 理 的 数据 进行 调 
归 模 型 ， 实 现 了 警觉 度 的 准确 估计 ， 并 结合 实例 数据 进行 了 整 ， 本 文 取 值 均 为 4。 
实验 验证 ， 证 明了 所 提 方 法 的 有 效 性 。 1.3 ”特征 降 维 
1 方法 通过 分 解 每 个 通道 的 EEG 信号 ,能够 提取 更 多 有 用 的 信 
g 息 ， 但 多 余 的 信息 也 可 能 造成 模型 训练 准 确 率 下 降 ， 曾 加 模 
1.1 特征 提取 型 训练 时 间 ， 因 此 对 脑 电 特征 数据 进行 降 维 处 理 有 利于 提高 
己 有 研究 表明 ，DE 在 人 员 警 觉 度 估计 中 的 性 能 优 于 传 警觉 度 估计 的 准确 率 ， 降 低估 计时 间 。 本 文 引 入 了 两 种 降 维 
统 功率 谱 密度 以 及 其 他 时 频 特 征 ， 因 此 本 文采 用 DE 作为 警 方法 : PCA 和 双向 二 维 主 成 分 分 析 (Two-directional two- 
觉 度 估计 的 特征 值 8204。DE 用 于 描述 连续 变量 的 复杂 性 ， 其 dimensionalPCA，TD-2DPCA)J05， 并 对 两 种 方法 在 警觉 度 估 
计算 公式 如 下 : 计 中 的 性 能 进行 了 对 比分 析 。 
DE=-| flog(f Cdr 0 传统 PCA 算法 是 一 种 利用 正 交 变 换 将 高 维 数据 映射 到 
、 低 维 空间 的 方法 ， 常 被 用 于 脑 电 数据 降 维 中 09， 能 够 在 一 定 
其 中 Fo 代表 随机 变量 的 概率 密度 函数 。 在 一 个 固定 的 频段 程度 上 降低 警觉 度 估计 时 间 , 但 PCA 在 对 高 维 数据 进行 降 维 
内 ，EEG 信号 的 分 布 可 以 近似 为 高 斯 分 布 Yo) ， 因 此 其 时 ， 需 要 先 将 高 维 数据 变换 为 行 向 量 ， 样 本 维 数 比 较 大 。 对 
DE 可 以 简化 为 于 NxMxL 维 的 脑 电 特征 ,在 每 个 时 间 窗 口 提取 的 EEG 特征 
2 0 是 一 个 NxM 维 的 矩阵 ， 不 同行 包含 的 是 不 同 脑 区 或 者 不 同 
De Be op Kier 0) 通道 的 信息 ,不 同 列 包含 的 是 不 同 子 频带 的 信息 , 而 PCA 在 
本 文采 用 两 种 DE 特征 提取 方法 : 一 是 对 于 每 一 个 通道 降 维 时 无 法 对 两 个 维度 的 信息 进行 区 分 。TD-2DPCA 能 够 直 
信号 , 均 在 5(0.5-4HZ)、9(4-8HZ)、Q (8-14HZ)、B(14-30HZ)、  ” 接 从 两 个 维度 对 原 有 和 矩阵 进行 降 维 ， 因 此 在 对 脑 电 的 通道 信 
7 (30-50HZ) 这 五 个 频段 上 提取 DE; 二 是 对 于 每 一 个 通道 信 息 和 频带 信息 进行 信息 提取 时 更 加 灵活 , 同时 TD-2DPCA 计 
号 ， 将 信号 的 某 一 部 分 频带 按照 一 定 的 频带 宽度 进行 等 分 ， 算 过 程 中 的 协 方差 矩阵 维度 远 低 于 PCA， 运 行 效率 更 高 ， 更 
比如 在 0-50HZ 的 频带 上 ， 按 照 2HZ 的 带宽 可 以 将 信号 分 解 适合 警觉 度 的 实时 估计 与 实际 应 用 
为 [0,2]、[2,4]、.…、[48,50] 共 25 个 频带 ， 然 后 从 各 个 子 频 带 在 2DPCA 中 ， 对 于 个 训练 样本 .1-12.,M ,每 个 样 
上 提取 DE。 所 有 DE 特征 都 是 采用 8 秒 长 无 重 铸 汉 宁 窗口 的 本 大 小 为 naxm 维 , 需要 将 其 通过 映射 矩阵 x 投影 到 主 成 分 空 
短 时 健 里 叶 变 换 来 计算 的 ， 因 此 得 到 的 特征 矩阵 维度 为 间 ， 即 : 
NxMxL，NN 代表 通道 数 ，M 代表 划分 的 子 频 带 的 数目 ，L Y=AX i=1,2,...,M (6) 
代表 整个 实验 中 提取 特征 的 数量 , 数值 上 等 于 EEG 信号 的 长 其 中 世代 表 第 ;个 样本 投影 的 结果 。 投 影 矩 阵 X 的 选取 准则 
度 除 以 窗口 的 长 度 。 是 使 得 投影 后 的 结果 辨识 度 更 高 ， 即 散 度 最 大 ， 因 此 定义 式 


1.2 ”特征 平滑 


在 警觉 度 估计 领域 , 大 部 分 研究 更 倾向 于 将 EEG 信号 映 F(X)=1(S.)=1{X" EA— EA) (A— FEA)X} 0) 
射 到 静态 离散 疲劳 状态 ， 但 实际 上 人 的 警觉 度 变化 是 一 个 时 其 中 : 5, 代表 Y 的 协 方差 矩阵 ， x(S 代表 协 方差 矩阵 的 迹 。 
间 动 态 过 程 四 ， 因 此 对 提取 的 特征 进行 平滑 不 仅 能 够 避免 测 ”定义 G = 可 4- 友 74- 友 为 样本 的 总 体 散 度 和 矩阵 , 可 以 发 现 
量 误差 带 来 的 影响 , 而 且 能 够 提高 特征 随时 间 变 化 的 连续 性 ， G, 是 一 个 维度 为 xxm 非 负 定 和 矩阵 ， 对 于 上 述 w 个 训练 样本 ， 
从 而 提高 人 员 警 觉 度 估 计 的 准确 率 。 滑 动 平均 是 常见 的 用 于 ”6 的 计算 方法 如 下 : 
数据 平滑 的 方法 ， 假 设 窗 口 长 度 为 ! 并 且 为 奇数 ， 对 于 脑 电 a (8) 
的 特征 序列 X ， 其 前 后 特征 的 权重 是 相等 的 ， 因 此 有 : 汉 页 

1 ， 1 志 < 
De i=(1+D)/2,...,L-(1-D/2 @) G -4 ND CN 一人) (9) 
五 else 计算 6 的 标准 特征 向 量 并 按照 特征 值 降序 排列 ， 由 于 
其 中 二 代表 原始 数据 ，z 代表 平滑 过 后 的 数据 ，L 为 待 处 理 。 是 一 个 维度 为 mxm 非 负 定 矩阵 ， 其 特征 向 量 有 mw 个 ， 则 投影 
数据 长 度 。 从 式 (3) 中 可 以 发 现 , 滑动 平均 无 法 对 [04+D/2-] 和 矩阵 的 最 优 解 为 X=[Xi,X,…,Xs] ， 其 中 X; 为 6 的 第 i 个 标准 
和 区 -4&-D/2+1 刀 区 间 的 原始 数据 进行 处 理 ， 这 可 能 导致 平 ”特征 向 量 。 投 影 结 果 了 =4X 的 维度 为 xxd ， 其 中 4 的 值 可 以 
滑 后 数据 在 两 端 出 现 跳 变 的 情况 ， 会 对 警觉 度 估 计 模 型 产生 ”通过 需要 的 信息 阔 值 设置 。 样 本 维度 从 ax 普 降 维 到 maxd ， 因 
较 大 的 影响 。 人 的 警觉 度 变化 是 一 个 随时 间 变 化 的 动态 特性 ， ”此 2DPCA 对 和 矩阵 进行 了 行 方向 的 投影 。TD-2DPCA 在 
在 短期 内 可 以 进行 估计 ， 基 于 此 对 传统 滑动 平均 方法 进行 了 2DPCA 的 基础 上 , 还 对 和 矩阵 进行 了 列 方向 的 投影 , 列 方向 的 
改进 : 投影 矩阵 z 求解 过 程 与 X 类 似 ， 只 需要 将 样本 4 进行 转 置 ， 

第 一 步 : 采用 传统 滑动 平均 对 数据 进行 平滑 ， 获 得 x ; 其 余 计算 过 程 相同 ， 最 终 得 到 的 样本 总 体 散 度 矩 阵 为 

第 三 步 : 从 zeoo 开 始 由 后 往 前 采用 变 窗 口 计 算 平 滑 结果 ， G, = E[(A— EA)(A— EA)’] (10) 


(7) 作 为 投影 


矩阵 Xx 的 目标 函数 : 
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售 度 为 ax 非 负 定 抢 阵 ， 


6, 是 


个 旨 
et rod pe ds 


Es 
C=2Z7AX (11) 
此 时 ，C 的 维 数 为 4xd ， 与 原 和 矩阵 维 数 (nxm) 相 比 ， 数 
据 得 到 有 效 压缩 ， 并 且 计 算 过 程 中 的 协 方差 矩阵 维度 分 别 为 


nxn 和 mxm 为， 相 比 PCA 计算 过 程 中 的 协 方差 矩阵 维 数 


(nm)x(nm) ， 计 算 效 率 得 到 提高 。 常 


于 控制 投影 结果 维度 的 


方法 有 两 种 : 一 是 直接 设 定 投影 后 的 


向 量 的 个 数 ， 二 是 根据 特征 值 计算 累 
率 第 一 次 大 于 设 及 定 阔 值 时 对 应 的 特征 
数 。 固 定投 影 结 果 的 维 数 会 影响 方法 


计 贡 献 率 可 能 出 现 某 一 特征 值 过 大 导 


维度 ， 即 设 定 选择 特征 
积 贡 献 率 ， 当 累积 贡献 
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其 特征 向 量 有 7 个 ， 驾驶 系统 中 进行 , 被 试 在 一 个 真实 的 车 辆 (没有 不 必要 的 引擎 

值 降序 排列 ， 得 到 列 方 向 的 投 和 其 他 组 件 ) 中 进行 长 时 间 驾 驶 , 车 辆 前 方 个 同步 变换 轰 
2， 将 样本 4 在 行 方向 和 列 方向 驶 场景 的 显示 屏 ， 道 路 主要 是 单调 直道 ， 更 容易 引发 疲劳 。 


一 共有 23 名 被 试 (12 


名 女性 ， 年 


有 被 试 都 保证 身体 健康 。 


龄 23.3/1.4) 参 与 了 实验 ， 所 
大 部 分 实验 在 13:30 左右 进行 ， 更 


容易 引发 被 试 疲劳 ， 


采集 了 18 个 通道 的 脑 部 EEG 信号， 
*， 分 别 是 FT7、 
TP8、 CPl、 CPZ、 CP2,\ Pl\ PZ\ P2、 


20 系统 ， 如 图 


一 次 实验 持续 时 间 约 2 小 时 。 实 验 一 共 


2 所 示 


有 极 放置 方法 为 国际 10- 
FT8、T7、T8、TP7、 
PO3、POZ、PO4、0O1、 


OZ、02， 其 中 CPZ 为 参考 电极 ， 因 此 数据 集中 包含 了 17 个 


值 个 数 即 为 投影 后 的 维 
的 泛 化 能 力 ， 而 设置 累 
致 降 维 后 信息 缺失 ， 因 


此 本 文 将 两 种 方法 进行 组 合 ， 整 体 上 
法 的 泛 化 能 力 ， 当 降 维 数据 维度 低 于 
固定 维度 降 维 。 
1.4 最 小 二 乘 支 持 向 量 机 


通道 的 EEG 信和 号 。 


5-7 号 电极 放 


位 置 为 SEED-VIG 数据 集 特有 ,如 图 3 


此 外 还 记录 了 7 个 通道 的 眼 部 和 
号 , 其 中 1-4 号 电极 放置 位 置 与 传统 测量 眼 电信 号 方式 


前 额 信 
样 ， 


| 


~ 
7 


了 请 “ 


该 信号 用 于 计算 


采用 累计 贡献 率 保证 方 


道 采集 的 信号 也 


最 小 允许 维度 时 ， 采 用 


SVM 在 小 样本 数据 的 分 类 和 预测 中 


有 较 快 的 运算 速 


度 和 精度 , 适合 警觉 度 的 实时 估计 ,LS-SVM 相对 SVM 计算 
更 加 简单 ， 运 行 速 度 更 快 ， 因 此 本 文采 用 LS-SVM ee 


特征 与 警觉 度 之 间 的 回归 模型 。 其 中 LS-SVM 的 核 函 数 采 
高 斯 核 函数 : 
Kn expC Ea) (12) 
其 中 ，o 代表 核 函 数 带 宽 。 
1.5 警觉 度 估 计 方 法 
警觉 度 估计 方法 的 核心 步骤 如 下 : 
Step1: 对 NN 个 通道 的 EEG 信号 ， 每 个 通道 信号 按 固定 


带宽 分 解 为 M 个 分 量 ; 
Step2: 
宁 窗 的 短 时 傅 里 叶 变 换 ，; 


在 每 个 分 量 上 提取 微分 焙 ,提取 方法 为 无 重 辣 汉 


Step3: 采用 改进 滑动 平均 和 TD-2DPCA 对 提取 的 微分 


炉 进行 平滑 处 理 和 降 维 ; 


Step4: 将 处 理 后 的 特征 输入 LS-SVM 回归 机 ,获得 人 员 


警觉 度 估计 结果 。 
本 文 方法 的 流程 图 如 图 1 所 示 。 
原始 EEG 固定 带宽 分 解 特征 提取 
We 
ON 
= 2 | 微分 炳 特征 
. re . 
Le 
v ¥ 
训练 集 测试 集 
1 1 
改进 滑动 平均 改进 滑动 平均 
y 了 
TD-2DPCA TD-2DPCA 
了 了 
警觉 度 标注 处 理 后 的 特征 处 理 后 的 特征 
汐 
-LS-SVM | 回归 模型 
了 
警觉 度 估计 结果 
图 1 警觉 度 估计 方法 框图 
Fig. 1 The framework ofthe vigilance estimation method 
2 ”结果 与 分 析 
2.1 脑 电 数 据 集 


数据 集 来 自 上 海 交 通 大 学 仿 脑 计 
公开 发 表 的 驾驶 


劳 数 据 集 SEED-VIG00。 


算 与 机 器 智能 研究 中 心 
实验 在 一 个 虚拟 


采样 到 200HZ， 
和 前 额 


言 号 降 采 档 


图 2 


眼 电 信号 (Electrooculogram, EOG), 4- 
可 以 视 为 EOG 污染 较为 严重 的 EEG 
为 了 降低 噪声 污染 ,减少 运算 量 ，17 个 通道 的 EEG 
并 进行 1-75HZ 的 带 通 
f 到 125HZ。 


国际 10-20 系统 电极 放置 图 


mm ar AN 
di va < 


重 滤 波 ; 7 个 通道 


es 
IE 
| 


Fpl FPZ FP2 


Fig.2 Electrode placement diagram of the International 10-20 system 


图 3 SEED-VIG 数据 集中 前 额 


2 


/ 


电极 放置 位 置 


Fig.3 Placement ofthe forehead electrodes in the SEED-VIG data set 


实验 采 / 


眼球 追踪 眼镜 记录 了 被 试 的 皮 眼 和 闭 眼 时 间 ， 


利用 PERCLOS 指数 进行 了 被 试 警觉 度 标注 。PERCLOS 指 


三 网 5 


为 8 秒 ， 
内 闭 眼 时 间 。 


PERCLOS = 


blink + CLOS 
time 


(13) 


中 ;time 代表 一 段 时 间 间 隔 , SEED-VIG 中 计算 指数 的 me 
blink 代表 时 间 间 隔 内 上 甩 眼 时 间 ，CZLOS 代表 时 间 间 隔 
因此 PERCLOS 指数 代表 了 被 试 在 一 段 时 间 内 
用 眼 和 闭 眼 时 间 的 占 比 ， 


取 值 范围 为 [0,1], PERCLOS 指数 越 


大 ， 代 表 被 试 警 觉 度 越 低 ， 疲 劳 程度 越 高 。 


2.2 ”基于 不 同 数 据 处 理 方 


法 的 警觉 度 估计 


本 文中 使 用 


到 的 信号 均 为 17 个 通道 的 EEG 信 


号 加 上 4 


个 通道 的 前 额 人 


十 
百 号， 


{ 21 个 通道 的 信号 ， 实 验 中 均 在 23 


个 被 试 数据 进行 了 实验 验证 。 


实验 中 对 


让， 获得 50 个 子 频 总 


种 平滑 方法 和 两 种 降 维 
比分 析 。 在 0-50HZ 频段 ， 以 1HZ 的 频带 宽度 对 信 
,然后 利用 短 时 傅 里 叶 变 换 在 每 个 频带 


担 


短 时 健 旦 


进行 特征 
的 DE 特征 维 


请 


多， 共 


的 性 能 进行 了 对 
号 进行 分 


方法 


有 E 叶 变换 中 的 窗口 为 8 秒 长 的 无 重 闭 汉 
宁 窗 。 数 据 集中 记录 了 118 分 钟 的 数据 ， 按 照 上 述 8 秒 长 窗 
了 885 个 数据 ， 
度 为 21x50x885, 代表 了 21 个 通道 , 50 个 频带 


因此 在 一 个 被 试 上 获得 
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和 885 个 样本 数量 。 提 取 的 特征 分 别 采 用 滑动 平均 、 改 进 滑 。” 每 次 实验 均 在 23 个 被 试 数据 上 了 进行 测试 。 
动 平均 两 种 方法 进行 平滑 处 理 ， 由 于 警觉 度 水 平 随 周期 长 度 实验 结果 如 表 2 所 示 。 从 Pearson 相关 系数 来 看 ，1HZ 
全 大 于 四 分 钟 H4， 因 此 滑动 平均 的 窗口 宽度 为 4 分 ”带宽 具有 最 高 的 精度 , 平均 数值 达到 了 87.45%, 从 RMSE 来 
钟 。 然 后 采用 两 种 降 维 方法 对 平滑 后 的 特征 进行 降 维 ， 为 了 看 ，2HZ 带宽 具有 最 高 的 精度 ， 平 均 数 值 在 9.07%， 表 明 与 
保证 实验 结果 具有 较 高 的 准确 率 同 时 数据 被 充分 降 维 ， 经 过 ”基于 3-5HZ 带宽 的 DE 提取 方法 相 比 ， 基 于 1-2HZ 带宽 的 
多 次 实验 实验 ， 将 PCA 的 贡献 率 设置 为 0.9，TD-2DPCA 在 ” DE 提取 方法 在 不 同 总 频带 上 都 具有 明显 的 优势 。Pearson 最 
行 方向 和 列 方向 的 贡献 率 分 别 设 置 为 0.9 和 0.85， 两 种 方法 ”大 对 应 的 数据 为 警觉 度 估 计 结 果 的 Pearson 相关 系数 最 大 时 
的 最 小 允许 维度 设置 为 4。 最 后 采用 LS-SVM 建立 特征 向 量 对 应 的 RMSE 和 总 频带 ， RMSE 最 小 对 应 的 数据 为 警觉 度 估 


起 


站 洒 


五 


与 PERCLOS 指数 之 间 的 回归 模型 ,其 中 x 和 o 设置 为 50000 计 结 果 的 RMSE 最 小 时 对 应 的 Pearson 系数 和 总 频带 ， 纵 向 
和 3000, 采用 均 方 根 误差 和 Pearson 相关 系数 估计 实验 结果 : 对 比 可 以 发 现 , 采用 1-2HZ 带宽 进行 警觉 度 估计 效果 优 于 其 
Ee 他 带宽 ， 并 且 基 于 1HZ 带宽 的 估计 结果 具有 最 高 的 Pearson 

WB (14) 。 精度， 为 91.50%， 基 于 2HZ 带宽 的 实验 结果 具有 最 优 的 


i RMSE 精度 ,为 8. ee Pearson 最 大 或 RMSE 
2 (15) 。 最 小 时 对 应 的 总 频带 范围 ， 发 现在 五 种 带宽 中 ， 最 高 估计 精度 
Fy- 二 (DJ - 二 (9) 村 的 总 和 几 大 多 有 和 是 从 OH 开始 ，20-30HZ 结束 ， 这 
表明 0-30HZ 的 EEG 中 包含 了 大 部 分 与 警觉 度 关联 的 信息 ， 同 

其 中 : yy 代表 测量 的 PERCLOS 指数 ，5 代 表 估 计 值 。N 代 ”时 进一步 表明 多 余 的 EEG 特征 提取 会 降低 警觉 度 估计 精度 。 
表 信 号 采样 点 数 。 实 验 中 采用 五 折 交 又 验 证 的 方法 ， 即 对 于 表 3 展示 了 所 有 实验 中 ， 估 计 结 果 的 Pearson 相关 系数 
每 一 个 被 试 ， 将 885 个 数据 分 为 五 段 ， 每 次 采用 其 中 四 段 作 大 于 0.9 时 对 应 的 RMSE， 以 及 特征 提取 中 采用 的 带宽 和 总 
为 训练 集 ， 剩 下 一 段 作 为 验证 集 ， 然 后 将 五 次 估计 结果 按 原 ” 频带 。 从 表 中 可 以 看 出 ， 一 共有 9 种 DE 提取 方法 能 够 实现 
有 顺序 排列 ， 得 到 六 ， 并 计算 RMSE 和 Pearson 相关 系数 。 较 高 精度 的 警觉 度 评估 ，RMSE 大 部 分 都 处 于 0.09 附近 。 同 
数据 平滑 和 降 维 对 警惕 度 估 计 的 影响 如 表 1 所 示 ， 表 中 时 ， 基 于 1HZ 带宽 的 方法 有 四 种 ， 基 于 2HZ 带宽 的 方法 有 4 
Pearson 相关 系数 和 RMSE 是 采用 23 组 被 试 实验 后 的 平均 中 ， 基 于 3HZ 带宽 的 方法 有 1 中 ， 基 于 4HZ 和 5HZ 带宽 的 方 
的 时 间 总 和 。 从 表 中 可 以 ” ”法 无 法 实现 较 高 精度 的 警觉 度 订 


Pearson, = 


值 , 时 间 是 23 组 被 试 测试 完毕 所 用 的 估 。9 种 DE 提取 方法 中 , 总 频 
发 现 ， 如 果 直 接 利用 原始 特征 向 量 进行 警觉 度 估计 ， 预 测 结 带 的 范围 没有 出 现 大 于 35HZ 的 情况 ， 这 表明 采用 DE 进行 警 
Pearson 相关 系数 只 有 64.08%， 并 且 运 行 的 时 间 最 长 ， 觉 度 佑 计时， 提取 DE 的 最 佳 总 带宽 在 0-35HZ 以 内 的 范围 。 
能 满足 实际 要 求 , 而 数据 经 过 平滑 和 降 维 后 ,， Pearson 相关 表 2 所 有 实验 结果 的 对 比分 析 
2 高 到 了 85% 以 上 ， 时 间 也 缩短 了 一 半 左 右 ， 因 此 对 Tab.2 Comparative analysis of all experimental results 
特征 进行 平滑 和 降 维 有 利于 警觉 度 的 实时 估计 ; 对 比 传统 滑 要 二 六 Pearson 最 大 RMSE 最 小 
动 平均 和 改进 滑动 平均 的 实验 结果 发 现 ， 改 进 滑动 平均 方法 轩 a Pearson RMSE 总 频带 RMSEPearson 总 频带 
性 能 优 于 传统 滑动 平均 方法 , Pearson 相关 系数 提高 了 约 2%， OWSE09 0 的 ED 0 (0 (HD 
RMSE 下 降 了 约 0.5%; 数据 经 过 改进 滑动 平均 和 TD-2DPCA 1 87.45 10.61 91.50 9.07 0-20 9.07 91.50 0-20 
处 理 后 进行 测试 ,得 到 的 Pearson 相关 系数 最 高 ,为 87.83%， 9 86.9 10.45 91.11 8.75 0-30 8.67 90.75 0-35 
运行 时 间 也 最 短 ，RMSE 比 最 低 值 仅 高 0.04%， 因 此 相对 于 3 86.26 10.6 90.41 922 0-25 9.1 89.87 4-20 
PCA 降 维 的 方法 ，TD-2DPCA 具有 一 定 的 优势 ， 能 够 应 用 于 4 85.77 11.03 89.91 9.77 0-25 9.08 89.76 0-15 
脑 电 特征 降 维 中 ， 并 且 随 着 特征 维度 的 增加 ，TD-2DPCA 降 5 85.41 10.9 88.25 9.65 425 9.15 87.73 4-20 
维 方 法 的 优势 会 更 加 明显 。 表 3 提出 方法 的 估计 结果 
表 1 不 同 数据 处 理 方法 及 结果 Tab. 3 The estimation results ofthe proposed method 
Tab. 1 Different data processing methods and results 带宽 (HZ) 总 频带 (HZ) Pearson(%) RMSE(%) 
处 理 方 法 Pearson(%) RMSE(%) 时 间 /s 1 0-20 91.50 9.07 
未 处 理 64.08 15.31 185 1 4-20 90.19 9.11 
滑动 平均 与 PCA 85.32 10.90 96 1 0-15 90.87 9.35 
滑动 平均 与 TD-2DPCA 86.02 11.21 94 1 0-25 90.01 10.09 
改进 滑动 平均 与 PCA 87.55 10.59 96 2 0-35 90.75 8.67 
改进 滑动 平均 与 TD-2DPCA 87.83 10.63 94 2 0-30 91.11 8.75 
2.3 基于 不 同 DE 提取 方法 的 警觉 度 估计 > 0-25 90.71 9.25 
为 了 进一步 提高 警觉 度 估计 精度 ， 本 文 对 DE 的 提取 方 2 0-20 90.13 9.28 
法 进行 了 详细 的 探究 。 分 别 将 总 频带 0-15 HZ、0-20 HZ、0- 3 0-25 90.41 9.22 
25 HZ、0-30 HZ、0-35 HZ、0-40 HZ、0-43 HZ、0-50 HZ、4- 2.4 与 现 有 方法 对 比 结果 
20HZ、 4-25 HZ、 4-30 HZ、 4-35 HZ、4-40 HZ、4-45 HZ、4- 为 了 验证 方法 的 有 效 性 ， 与 现 有 的 警觉 度 估计 方法 进行 


50HZ、 8-25 HZ、8-30 HZ、8-35 HZ、8-40HZ、8-45 HZ、8- 了 对 比 , 结果 如 表 4 所 示 。 其 中 提出 方法 从 2HZ 带宽 分 解 0- 
50 HZ， 分 别 以 1HZ、2HZ、3HZ、4HZ、5HZ 的 带宽 进行 等 35HZ 总 频带 后 的 分 量 提取 微分 焙 ， 五 段 带 方法 从 五 频带 中 
分 ， 然 后 从 获得 的 子 频带 上 提取 DE。 如 将 0-25HZ 频段 以 提取 微分 炉 ， 文 献 [12] 是 从 2HZ 带宽 分 解 0-50HZ 总 带宽 后 
2HZ 进行 等 分 ， 可 以 获得 0-2HZ、2-4 HZ， 直 至 22-24 HZ 共 的 分 量 提取 微分 粹 。 从 表 4 可 以 看 出 ， 提 出 方法 的 警觉 度 估 
12 个 子 频带 , 其 他 频段 的 等 分 方法 与 之 类 似 。 最 终 一 共有 105 ， 计 结果 为 0.91/0.09, 文献 中 警觉 度 估计 精度 最 高 为 0.85/0.09， 
种 DE 提取 方式 。 特 征 提 取 采 用 短 时 传 里 叶 变 换 ， 数 据 平滑 ”表明 提出 方法 的 性 能 优 于 现 有 方法 。 同 时 对 比 提出 方法 和 五 
采用 改进 滑动 平均 方法 ,数据 降 维 采 用 TD-2DPCA 方法 , 下 段 带 方 法 的 估计 结果 可 以 发 现 ， 在 特征 提取 过 程 中 ,将 EEG 
归 机 采用 LS-SVM, 所 有 算法 的 参数 设置 与 第 一 个 实验 相同 ， 细 分 后 提取 特征 比 从 五 频带 提取 特征 更 有 效 。 
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表 4 与 现 有 方法 对 比 结果 


Tab.4 The comparison results with existing methods 


方法 特征 处 理 归 器 Pearson RMSE 


日 


提出 方法 改进 滑动 平均 +TD-2DPCA SVM 0.91 0.09 


五 频带 方法 “改进 滑动 平均 HTD-2DPCA SVM 0.89 0.09 
滑动 平均 SVM 0.83 0.1 
Zheng 等 滑动 平均 CCRF 0.84 0.1 
滑动 平均 CCNF 0.85 0.09 


3 


对 


了 系统 的 研究 ， 提 出 了 一 种 基于 改进 滑动 平均 、TD-2DPCA 


结束 语 


为 提高 基于 脑 电 DE 特征 的 警觉 度 实时 估计 精度 ， 本 文 
DE 的 提取 方法 ， 数 据 的 平滑 方法 和 数据 的 降 维 方法 进行 


和 最 


EDE 提取 的 警觉 度 估计 方法 。 首 先 结合 警觉 度 的 时 间 


动态 特性 ， 利 用 变 窗口 平均 和 短 时 预测 的 方法 对 传统 滑动 平 


均 方 法 进行 了 改进 ， 提 高 了 边缘 数据 与 整体 数据 的 关联 ; 在 
数据 降 维 上 ， 引 入 TD-2DPCA， 对 比 传统 PCA 方法 ,运算 速 
度 和 信息 提取 能 力 都 有 小 幅度 的 提高 ， 最 后 结合 实例 数据 和 
改进 的 方法 ， 详 细 对 比分 析 了 不 同 总 频带 和 不 同 分 解 带 宽 对 
警觉 度 估计 的 影响 ,结果 发 现 ,总 频带 在 0-35HZ 以 内 , 分 解 
带宽 在 1HZ 或 2HZ 时 ， 警 觉 度 估计 的 精度 较 高 ，Pearson 相 
关系 数 约 为 0.91，RMSE 约 为 0.09， 并 且 优 于 现 有 方法 的 估 
计 精 度 0.85/0.09。 下 一 步 将 结合 迁移 学 习 算 法 ， 将 提出 方法 


MY 
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